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ABSTRACT 

We present deep WIYN Ha imaging of the dwarf irregular starburst galaxy NGC 
1569, together with WIYN SparsePak spatially-resolved optical spectroscopy of the 
galactic outflow. This leads on from our previous detailed analyses of the state of 
the ISM in the central regions of this galaxy. Our deep imaging reveals previously 
undetected ionized filaments in the outer halo. Through combining these results with 
our spectroscopy we have been able to re-define the spatial extent of the previously 
catalogued superbubbles, and derive estimates for their expansion velocities, which we 
find to be in the range 50-100 kms^^. The implied dynamical ages of <25 Myr are 
consistent with the recent star- and cluster-formation histories of the galaxy. Detailed 
decomposition of the multi-component Ha line has shown that within a distinct re- 
gion ~700 X 500 pc in size, roughly centred on the bright super star cluster A, the 
profile is composed of a bright, narrow (FWHM < 70 kms~^) feature with an under- 
lying, broad component (FWHM ^ 150 kms^^). Applying the conclusions found in 
our previous work regarding the mechanism through which the broad component is 
produced, we associate the faint, broad emission with the interaction of the hot, fast- 
fiowing winds from the young star clusters with cool clumps of ISM material. This 
interaction generates turbulent mixing layers on the surface of the clouds and the 
evaporation and/or ablation of material into the outflow. Under this interpretation, 
the extent of the broad component region may indicate an important transition point 
in the outflow, where ordered expansion begins to dominate over turbulent motion. 
In this context, we present a multi- wavelength discussion of the evolutionary state of 
the outflow. 

Key words: galaxies: individual (NGC 1569) - galaxies: starburst - galaxies: ISM - 
ISM: kinematics and dynamics - ISM: jets and outflows. 



1 INTRODUCTION 

Outflows powered by the collective injection of kinetic en- 
ergy and momentum from massive stars and supernovae 
(SNe) can drastically affect the structure and subsequent 
evolution of galaxies. Thus, a good understanding of the 
feedback mechanisms between massive stars, star clusters 
and the ISM is fundamental. In particular, dwarf galaxies 
are thought to be strongly affected by the effects of feed- 
back since their smaller gravitational potentials mean that 
supernova- heated gas can escape more easily ( Larson] 1974 1. 
Although the ejection of the ISM potentially could have 
significant consequences on the star-formation rate within 



out may not be as eflicient as first thought (De Young & 
|Heckman||l994[ |Martin||1998[ ). It is therefore important to 
study such systems to understand how gas is removed and 
what effects this has in the evolution of the galaxy. 

NGC 1569 (UGC 3056, Arp 210, VII Zw 16, IRAS 



4260-1-6444) is a nearby (2.2 Mpc; |Israel||1988|l, low met- 
alicity (0 25 ZQ ; |Devost, Roy fc Drissen|1997||Kobulnicky"fc 
Skillman|1997t dwarf irregular galaxy that has recently un- 
dergone a period of starburst activity. The most recent burst 
is thought to have peaked between ~10-100 M yr ago with an 



average star-formation rate of ~0.5 Mq yr ( Greggio et al 



19981. At some point near the end of this event, the two 



these low-mass systems (Dekel & Silk 19861, more recent well-known super-star clusters (SSCs) A and B were formed 



work suggests that ejection of hot gas through bubble blow- 



( Arp & Sandage 


1985 


O'Connell, Gallagher & Hunter|1994 


de Marchi et al. 


1997 


1, and together with the slightly older 
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to dominate the energetics of the central regions. 
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H I observatio ns of NGC 1569 | |Stil fc Israel|1998|[2002 
Miihle et al.|2005 1 clearly show morphological and kinemat- 



ical signatures caused by the starburst. Firstly, the Hi kine- 
matics are "strongly disrupted" within the central 900 pc, 
with little or no evidence for the disc rotation seen at lager 
radii. Furthermore, large seemingly tidal structures are seen, 
including a so-called bridge connecting the galaxy to a low- 



mass Hi cloud (the companion; Stil &: Israel||1998 1, a large 
H I arm extending to the west of the disc, and a very faint 
filamentary Hi stream wrapping around the south of the 



disc (Miihle et al. 20051. A 'hot spot' (a region of veloc- 
ity crowding) on the western edge of the disk was found by 
Miihle et al. (2005), who interpreted it as the impact loca- 



tion of infalling gas from the companion. This provides a 
compelling explanation of how the starburst event was trig- 
gered. 

Ha images of this galaxy show an equally chaotic, com- 
plex structure to the warm ionized component, exhibiting 
many filaments, bubbles and loops. Many of these were iden- 



tified by Hunter, Hawley & Gallagher ( 1993 1 from deep Ha 
imaging. Later kinematical studies found that these fila- 
ments form part of a cellular outflow structure, compris- 
ing large-scale expanding superbubbles on the northern and 
southern sides of the disc ( Tomita, Ohta fc Saito|1994 Heck- 



man et al. 11995 Martin 19981. By analysing spectra from 



two long-slits placed parallel to the major and minor axes. 



Heckman ct al. ( 1995 1 found evidence of shocks in the outer 



regions of the ionized halo. The ratios of [Ol]/Ha, [Sll]/Ha 
and [Nil] /Ha (used often to diagnose and trace the condi- 
tions within ionized gas; |Veilleux fc Osterbrock|1987||Dopita| 
et al.|2000 l were all found to be high in the outer filaments. 



indicating either an increase in the importance of shocks, or 
that photoionization becomes less important with distance 
as the ionizing radiation from the central starburst becomes 
diluted. 

The existence of shocked gas is supported by high- 
resolution X-ray observations of NGC 1569. [Martin, Kob-| 
[ulnicky fc Heckman] ( |2002[ ) examined the X-ray properties 
of the outflow with Chandra, and found significant soft (0.3- 
0.7 keV), diffuse emission coincident with the Ha morphol- 
ogy. Although they found the X-ray colour variations to be 
inconsistent with a free-streaming wind, they concluded that 
the X-ray emission probably originates in the halo shock 
generated by the outflowing gas, possibly from the mixing 
layers between the shock and the bubble interior. 

In order to properly characterise the structure of the 
outflow, it is essential to study its kinematics. From opti- 



cal long-slit spectroscopy, Martin ( 1998 1 found the outflow 



to be composed of numerous superbubbles. In general, she 
finds the redshifted component of the split-line profile to be 
stronger in the south and weaker in the north, suggesting an 
inclined outflow aspect. This is consistent with more recent 
X-ray absorption measurements ( jMartin et al.|2002j ) and H I 
observations (from which an inclination angle of ~60° was 
derived [Stil fc Israel||200"2| ). Although it is unclear whether 
the two SSCs, A and B, alone are sufficiently powerful to 
provide enough mechanical energy and ionizing radiation to 
drive the whole outflow and power the galaxy's diffuse ion- 
ized medium ( Martin|1997 Martin et al.|200^ I , what is clear 
from the Ha morphology is that energy is being injected 
throughout the central starburst zone of the disc from mul- 
tiple sources. 



A detailed look at the spectral line proflles, however, 
shows that near SSC A, "distinctly non-Gaussian" Ha emis- 
sion can be found, exhibiting weak but very broad wings 
(Heckman et al. 1995| ). Broad emission line wings have 



been detected in other nearby starburst galaxies (e.g. 



Izo- 



tov et al.||1996| [Homeier fc G allagher| |1999| [Marlowe et al 
19951 [Mendez fc Esteban[|1997: Sidoli, Smith fc Crowther 



2006 , Westmoquett e et al.|2007 ,) . Due to mismatches in spec- 
tral and spatial resolution and in the speciflc environments 
observed, the nature of the energy source for these broad 
lines has been contested, and has resulted in the proposal of 
number of possible explanations. However, a detailed IFU 
(integral field unit) study of the ionized ISM conditions in 
four regions within the 200 pc region surrounding SSC A by 



Westmoquette et al. ( 2007 hereafter Paper I \ and Westmo- 



quette et al. ( 2007 hereafter Paper II I have shed a consid- 



erable amount of light on this problem. 



By accurately decomposing the emission line profiles 
across each of the IFU fields, we found the line shape to be, 
in general, composed of a bright narrow feature (intrinsic 
FWHM ~ 50 kms~^) superimposed on a fainter broad com- 
ponent (FWHM ~ 200-400 kms"^). By mapping out their 
individual properties, we identified a number of correlations 
between the line components that allowed us to investigate 
in detail the state of the ionized ISM. We concluded that 
the broad underlying component is most likely produced in 



a turbulent mixing layer (Slavin et al. 1993 Binette et al 



19991 on the surface of cool gas knots, set up by the impact 



of the fast-flowing cluster winds | |Pittard et al.||2005) . Our 
analysis revealed a very complex environment with many 
overlapping and superimposed components, but surprisingly 
no evidence for organised bulk motions. We concluded that 
the four regions sampled are all located well within the wind 
energy injection zone ( Shopbell fc Bland-Hawthorn] 1998 1 at 
the very roots of the outflow, and that the collimation pro- 
cesses required to transform the turbulent motions into an 
organised net outflow forming the large-scale superbubbles 
must occur between 100-200 pc from the central star clus- 
ters. 



With this in mind, we have obtained new deep Ha imag- 
ing of NGC 1569, together with spatially-resolved spectra of 
the outer halo, to investigate in detail the morphology and 
kinematics (including the line proflle shapes) of the warm 
ionized component of the outflow at these large radii. In 
Section [2] we present the observations and in Section [S] we 
map out the properties of the line components and discuss 
the results of the spectroscopy. Since a number of SparsePak 
fibres are coincident with or lie adjacent to the Gemini 
GMOS/IFU fields presented in Papers I and II, we compare 
the results obtained with the two instruments in Section [J] 
In Section|5]we discuss the state of the outflow, including the 
conditions in the inner and outer halo, and we summarise 
our flndings and conclusions in Section |6] 



The galactic outflow in NGC 1569 3 




N 



50" ~ 500 pc 



Figure 1. Continuum subtracted Ha image of NGC 1569 taken witli WIYN/MiniMo. Tlie intensity stretch has a cube-root scaling and 
is inverted to enhance low surface-brightness features. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Imaging 

We obtained WIYKQs.S-m Mini-Mosaic (MiniMo) imaging 
of NGC 1569 on 13th December 2006 through Harris J?-band 
and W15 (= Hq) filters, in seeing conditions of ~l'.'l-l'.'2. 
Mini-Mos consists of two 4096 x 2048 CCDs separated by a 
small gap. Two offset exposures were taken for each filter in 
order to cover the gap and to aid cosmic ray and defect re- 
moval; exposure times were 300 s for each i?-band image and 



^ The WIYN Observatory is a joint facility of the University of 
Wisconsin-Madison, Indiana University, Yale University, and the 
National Optical Astronomy Observatories. 



1200 s for each Ha image. Dome fiats and bias frames were 
also taken. Reduction was achieved using the NOAO MS- 
CRED package within iRAlj^ First, bias and fiat fields were 
combined and subtracted/divided into the object frames, 
then an accurate world coordinate system was applied us- 
ing the tasks msczero and mscmatch. Single images were 
constructed from the individual chips using mscimage, and 
large-scale gradients in the sky background were removed us- 



^ The Image Reduction and Analysis Facility (iraf) is dis- 
tributed by the National Optical Astronomy Observatories which 
is operated by the Association of Universities for Research in 
Astronomy, Inc. under cooperative agreement with the National 
Science Foundation. 
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50"- 500 pc 2) 

Figure 2. Hq morphology of the NGC 1569 halo. Numbered filaments correspond to the nomenclature of [Hunter et al.] l |19 93l|; newly 
identified filaments are numbered from 14 onwards. The arc 1/arc 2 nomenclature is fro m |WalIer| ( [1991^ . The dashed regions labelled 
with letters correspond to large-scale superbubblcs and follow the nomenclature of [Martin] l |1998[ ). The white cross marks the position of 
SSC A. 



ing MSCSKYSUB. Finally the separate mosaiced frames were 
combined using mscstack, after matching their overall in- 
tensities using MSCIMAGE. The final mosaiced R and Ha im- 
ages cover a field-of-view of 9.5 x 10.7 arcmins at a spatial 
sampling of 0'.'14 pixel" ^. 

The 7i-band image was used as a proxy ofT-band ex- 
posure to subtract the continuum from the Ha image. The 
final continuum-subtracted Ha mosaic is shown in Fig. [l] 
Although the depth of this image equals that of the deepest 
image previously published ( Hunter et al.|19"93 l, its quality 
(particularly around the bright saturated star in the north of 
the image) and dynamic range allow previously undetected 
features to be seen. For the first time, the morphological 
structure of the Ha-emitting gas can be easily traced from 
the outer halo right into the central regions, where multi- 
ple nested bubbles and shells are organised in a seemingly 
chaotic manner. Sub-structure can be seen down to the res- 
olution limit of the observations, and emission is discernible 
out to distances of 2.5-3 arcmins (1500-1800 pc) from the 
disc. The halo morphology seen in Fig.[l]is clearly very com- 
plex. 



2.1.1 Shell morphology 

Kinematical studies of the NGC 1569 halo have identified 
a number of large-scale supershells through their double- 
peaked signatures in bright emission lines. Through these 



methods, Martin ( 1998 1 found the southern flow to be com- 



posed of three main expanding shells (which she denoted 
F, G, and A). The largest line-splitting was found in shell 
A, outlined on the western side by the prominent 'western 
arm' ( jHodge|1974 1, corresponding to expansion velocities of 
Wcxp ~ 70-100 kms~^ and indicating that this bubble might 
not be "strongly decelerated as it pushes through galactic 
halo" . The northern halo was also found to be composed of 
three large-scale shells (E, D and B) where the expansion ve- 
locity of shell B was found to be Uexp ~ 80 kms"^. Shell C, 
also identified by Tomita et al. ( 1994 1 in their stepped long- 



slit observations, represents the central cavity in NGC 1569 
centred on SSCs A and B, and measures ~20" in diameter 
with a Uoxp of ~40 kms~^. 

Deep Ha imaging reveals a cornucopia of structure 
within the halo that was first catalogued by [Hunter et al.| 
( ]1993 1. The depth and quality of our new imaging (Fig. [l| 
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shows many more features that can now be added to the hst. 
In Fig. |2] we have labelled the prominent Hq halo features 
following the nomenclature of Hunter et al. (19931, where 
new detections are numbered from 14 onwards. To clarify 
matters, we suggest that: filament 4 refers to the whole set 
of linear streamers emanating from the north-west of the 
disc; filament 5 refers to the bubble-like curved filament, 
whereas filament 18 refers to the brighter linear streamer; 
and filament 6 refers to the whole of the western arm. 

Based on the imaging and spectroscopy presented in 
this paper, we have re-examined the existence and extent of 



the large-scale superbubbles catalogued by [Martin ( 1998 1 



Previously undetected faint structure and spectroscopically 
identified expanding material (see Section [sj have allowed 
us to refine the spatial extents of supershells A and B. How- 
ever, the Ha morphology does not support a distinction be- 
tween shells G and F (Martin|1998 1, so we propose that the 
whole south-eastern shell complex consisting of filaments 5, 
7, 8, 17, 18, 19, and 20 should be referred to as superbubble 
F. This bubble corresponds to the south-eastern X-ray lobe 
( Martin et al.|2d02 and Section 5.21. In our interpretation, 
the two kinematically identified expanding structures that 



led Martin (19981 to separately define shells G and F, rep- 
resent just one level of a 'hierarchy of structure', where mul- 
tiple nested bubbles (e.g. 19 and 8 or 17 and 5) are contin- 
uously interacting to form larger structures. This interpre- 
tation can also be applied to the north-eastern halo, where 
a distinction between shells D and E is equally difficult to 
corroborate from the Ha and X-ray morphologies (see Sec- 
tion 5.2 1. Thus we suggest superbubble E should encompass 
the whole north-eastern halo, including both the faint Ha 
wisps in this region (evident from both the Hq imaging and 
spectroscopy presented below) and the extent of the north- 
eastern X-ray lobe (Martin et al. 2002[ and Section 5.2 I. The 
spatial extent of the large-scale superbubbles according to 
our interpretation are drawn on Fig. [2] with dashed lines and 
are labelled with the letters A, B, E and F. 

The inclination of the disc (~60°; |Stil fc Israel|200"2| is 
such that the northern side is tipped away from our sight- 
line. This scenario is supported by the consistent reversal of 



line strengths between the north and south ( Heckman et al 



|1995t|Martin|1998 


\ and the more homogeneous H I distribu- 


tion in the north ( 


Miihle et al.||2005L Therefore, unlike the 



southern side where we can, in effect, see up into the disc, 
the northern side is much more obscured. A large number of 
linear streamers (labelled collectively as filament 4) all point 
away from the central region containing SSCs A and B, and 
appear to form the tails of the swept-back, cometary-like 
shape of the gas knots seen just to the north of SSCs A and 



B ( Paper II \ . These radial Ha filaments could be compressed 
regions of shells or ablation trails from dense gas clouds (e.g. 
Melioh et al.|2005 l. 



2.2 Spectroscopy 

We observed four regions surrounding NGC 1569 with 



the SparsePak instrument ( Bershady et al. 2004 1 on the 



WIYN telescope during the period 13-16th December 2004. 
SparsePak is a "formatted field unit" similar in design to 
a traditional IFU, except that its 82 fibers are arranged in 
a sparsely packed grid, with a small, nearly-integral core 



numbers). SparsePak was designed to maximise through- 
put and spectral resultion at the expense of spatial cover- 
age/resolution, and as such has a total light throughput of 
~90 per cent longwards of 500 nm. Each fibre has a diameter 
of 500 \im, corresponding to 4'.'69 on the sky; the formatted 
field has approximate dimensions of 72 x 71.3 arcsecs, in- 
cluding seven sky fibres located on the north and west side 
of the main array separated by ~25". The mapping order 
of fibers between telescope and spectrograph focal planes 
was purposefully designed in a fairly randomised fashion in 
order to distribute sky fibres evenly along the slit and min- 
imise the effects of spectrograph vignetting on the summed 
sky spectrum. SparsePak is connected to the Hydra bench- 
mounted echelle spectrograph which uses a 2048 x 2048 CCD 
detector. Using an order 8 grating with an angle of 63.25° 
(giving a spectral coverage of 6450-6865 A and dispersion 
of 0.20 Apix~^), we were able to cover the nebular emis- 
sion lines of Ha, [N ii] AA6548, 6583, and [S ii] AA6716, 6731. 
A number of bias frames, flat-fields and arc calibration ex- 
posures were also taken together with the science frames. 

The positions of the SparsePak fibres for all four fields 
are shown in Fig. |3j overlaid on the WIYN Ha image. For 
clarity, we have named each position with a number, and a 
bracketed letter referring to the name of the supershell with 
which the position is associated. The arrangement of the 
fibres (including the offset sky fibres) can clearly be seen. 
A list of the coordinates for each position, and the total 
exposure times are given in Table [l] The four pointings were 
chosen to sample the outer galactic wind flow, and coincide 
with the main expanding shells identifled by [Heckman et aT] 
( 1995 1 and Martin ( 1998 1 and with regions of bright diffuse 
X-ray emission ( [Martin et al.|2002[ ). 

Basic reduction was achieved using the NOAO hydra 
package within iraf. The first step was to run APFIND on 
the fiat-field exposure to automatically detect the individ- 
ual spectra on the CCD frame. The output of this task is an 
aperture identification table which can be used for each sci- 
ence frame to extract out the individual spectra. The task 
DOHYDRA was then used to perform the bias subtraction, 
flat-fielding and wavelength calibration. The datafiles at this 
stage contained 82 reduced and wavelength calibrated spec- 
tra, each 1 pixel in width and ordered by the position they 
were recorded on the CCD (the 'fibre order'). Because of 
the positioning of our fields on NGC 1569, many of the sky 
fibres fell on regions of strong nebular emission (see Fig. [s]) 
and were therefore unusable. By a combination of visual in- 
spection of both the finding chart and the sky spectra them- 
selves, we determined which fibres to use for the sky sub- 
traction process (those that were contaminated were used 
as normal science fibres). A list of which fibres were used 
for sky subtraction in each position is given in Table [l] Af- 
ter determination of the appropriate usable sky fibres, an 
average sky spectrum was created by averaging the spectra 
and sky-subtraction was achieved by re-running dohydra 
with the sky-subtraction option switched on. Cosmic-rays 



were cleaned from the data using lacosmic (van Dokkum 



(see Bershady et al.[[2004 their fig. 2 for a map of the fibre 



2001 1, before final combination of the individual frames was 
achieved using imcombine. The final datafiles now each con- 
tained 82 reduced, wavelength calibrated and sky-subtracted 
spectra. 

In order to determine an accurate measurement of the 
instrumental contribution to the spectrum broadening, we 
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Figure 3. WIYN Ho image of NGC 1569 with the positions of the four SparsePak fields (including sky fibres) overlaid. The letters in 
brackets indicate the main supershell that they cover. The bold lines indicate the extent of the broad component region as described in 
Section [5.1| The white cross marks the position of SSC A. Inset: HST /'WFPC2 F656N image covering the white dashed area outlined on 
the main image, showing the position of the Gemini GMOS/IFU fields of view (5 X 3.5 arcsecs; [Paper I| [Paper on an HST /WFPC2 
F656N image, compared to the nearest SparsePak fibres. 



Table 1. WIYN SparsePak observations 



Position 


RA Dec 


Exposure Time 


Sky fibres used 




(J2000) 


(s) 


for subtraction 


1(E) 


4h3imoi!0 64°51'25yo 


8 X 1800 


16,22 


2(B) 


4'^30™48!5 64°51'3r'0 


9 X 1800 


2,37,54,80 


3(F) 


4h30™52?0 64°49'18^'0 


8 X 1800 


70,80 


4(A) 


4'^30™42!0 64°50'15yo 


8 X 1800 


2,37,54,80 



selected high S /N spectral lines from a wavelength calibrated 
arc exposure that were close to the Ha line in wavelength, 
and sufficiently isolated to avoid blends. After fitting these 
lines with Gaussians in all 82 apertures, we find the average 
instrumental width is 30.9 ± 0.4 kms~^. 



2.2.1 Line fitting 

We fitted model Gaussian profiles to each of the emission 
lines detected in each spaxel of the four Spars ePak fields us - 
ing the fitting routine PAN (Peak ANalysis; Dimeo|2005 1. 
A detailed description of the program and the customisa- 
tions we have made to it are given in [PaperT] In order to 
maintain a consistent approach to help in the minimisation 
of the fit and in the analysis, the first component of the fit 
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Figure 4. Example Ha line profiles chosen to represent the main types of profile shapes observed over the four fields (top four) and to 
aid explanation in the text (lower two), (a) split line profile; (b) narrow, bright component with underlying, broad component; (c) split 
narrow profiles with underlying broad component [from N-E of position 4(A)]; (d) triple component line [from centre of shell B, position 
2(B)]; (e) a spaxel from near the south-east of position 2(B) (30" from SSC A) that does not show a broad underlying component, only 
two narrow components; (f) a low S/N line broadened by multiple unresolved components (from the north-west of position 2(B) on the 
edge of shell B). Observed data is shown by a solid black line, individual Gaussian fits by dashed red lines (including the straight-line 
continuum level fit), and the summed model profile in solid red. Flux units are arbitrary but on the same scale. 



initial guess (hereafter referred to as CI) was always assigned 
to the brightest of any multiple-components encountered. In 
many cases, a second Gaussian component (hereafter C2) 
was needed; sometimes it was required to fit a broad under- 
lying profile, and in others it was required to fit the fainter 
component of a split line. The number of components fitted 
to each line was determined using a combination of visual 
inspection and the fit statistic output by PAN. If, as was 
often found in the outer parts of the galaxy where the S/N 
of our spectra were lower, only one component was detected, 
but its velocity was consistent with the surrounding C2 ve- 
locities, we manually re-assigned the component to C2 in 
order to aid interpretation. 

A range of example Ha line profiles from various spaxels 
are shown in Fig. |4j together with the individual Gaussian 
fits required to model the integrated line-shapes. The top 
four profiles were chosen to illustrate the main types of Ha 
line shapes encountered over the four fields; the majority are 
similar to the top two examples, being composed of single- or 



double-Gaussians (most being simply double-peaked - panel 
a - but a good proportion exhibiting a broad underlying 
component - panel b). In only five fibres do we find a clear 
triple profile in Ha, two of these being of the type shown 
in Fig. |4|; (three distinct narrow components), and three 
being of the type shown in panel d (broad underlying and 
faint secondary narrow component). The bottom two panels 
(e and f) show example profiles that are discussed in the 
following sections. 



3 EMISSION LINE MAPS 

The spatial distribution of the properties of each Gaussian 
component identified in the SparsePak fields were mapped 
using the visualisation tool, DAISY, described in [Paper I] 
Figs [5] [6] and [T] show respectively the flux, FWHM and ra- 
dial velocity maps of the Ha line for the four SparsePak po- 
sitions, overlaid on a faded reproduction of the WIYN Ha 
image. For clarity, we do not include the third component 
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since it is only detected in a few spaxels. In CI (the bright- 
est component) we detect Ha emission out to the extent of 
our spatial coverage - a maximum of 2.4 arcmins (1.5 kpc) 
away from SSC A. 

Fig. [5] shows the flux maps of CI (left panel) and C2 
(right panel) seen by SparsePak. It is clear that the flux vari- 
ations of both components correspond well to the Ha image 
shown in the background. The average flux ratio of C2/C1 
over the four fields is ~0.7 ± 0.4. This value gives an esti- 
mation of the importance of the different components; un- 
fortunately the quality of these datasets is not high enough 
to reveal any trends in this value across the observed ar- 
eas. Our derived C2/C1 ratio is considerably higher than 
that found from any of our CMOS IFU fields of the cen- 
tral regions (0.1-0.3; Paper Il| . The reasons for this will be 
discussed in Section |4] 

In Fig. [6j we plot the spatial variation of the FWHM 
of Ha. The line width of CI (left-hand panel) remains fairly 
constant across the face of the galaxy, ranging over ~20- 
70 kms~^ (corrected for the instrumental contribution). In 
a number of places broad CI lines are measured, but be- 
cause these are almost all located adjacent or near to spax- 
els where multiple components are detected, these broad fits 
are likely to be caused by the convolution of ~two unre- 
solved profiles. An example of a broad line from one of the 
outer-wind regions is shown in Fig. [4j, illustrating how the 
low S/N has meant that only a single-Gaussian is required 
for a satisfactory fit. However, the width of this Gaussian 
suggests that the profile contains multiple kinematic com- 
ponents. The FWHM map for C2 is shown in the right-hand 
panel of Fig. [6] At large radii, where a second component is 
detected it typically has a width similar to that of CI, and 
forms part of a classic split-line profile (e.g. Fig|4|i). The 
broadest lines (FWHM < 160 kms~^) are found within the 
central regions of the galaxy. Here C2 represents an under- 
lying broad component of the type seen in Fig. ^p, and very 
similar to what we found near the centre of the galaxy from 
our CMOS IFU study ( [Paper I| [Paper II| ). The variable S/N 
nature of the observations makes it difficult to determine 
whether at larger radii we stop detecting this broad compo- 
nent because it actually does not exist, or whether as the Ha 
emission becomes weaker with increasing galactocentric dis- 
tance, it fades into the continuum noise. Detailed inspection 
of the line profiles show that we can identify the broad com- 
ponent within a definable region (hereafter referred to as the 
'broad component region'), 70 x 50 arcsecs (~700 x 500 pc), 
approximately centred on SSC A and roughly aligned with 
the galaxy major axis. This region is outlined with a solid 
line in Fig. [3j 

The fact that this broad component region has an ex- 
tent approximately equivalent to the Ha-bright disc could 
suggest that the broad emission simply fades into the con- 
tinuum at this point and becomes undetectable. However we 
find evidence supporting a definite physical extent to the ex- 
istence of the broad component. Fig. [4^ shows an example 
Ha profile from just outside this broad component region 
(at a distance from SSC A of « 330 pc), and clearly shows 
how it can be modelled with two Gaussians of similar narrow 
widths, with no requirement for a broad underlying compo- 
nent. Correspondingly, we find spaxels in a number of places 
just within the border of the outlined region with an equiv- 
alent S/N as the example shown, but clearly exhibiting the 



Table 2. Measured and derived supershell properties 



Supershell 




Diameter 


Dynamical 




(kms"-') 


(kpc)" 


age (Myr)" 


A 


~90 




< 10-15 


B 


~85 




< 10-15 


E 


<50 


~1.3 


<25 


F 


<100 


~1.3 


<25 



" see Section Hm 

broad-line component. Therefore, we interpret these results 
to mean that the broad line does physically cease to exist 
at the boundary of this region. This finding has important 
consequences that we discuss in Section [5. 1| 

Fig. [7[ shows the spatial variation of the Ha radial ve- 
locities with respect to the systemic velocity of the galaxy, 
for which we have adopted the value of Usys = —80 kms~^ 



(Paper I and references therein). In general, the velocity of 
CI (left-hand panel) appears consistent with the systemic 
velocity, and remains fairly constant over the field, varying 
between —40 and 4-40 kms~^ relative to Vsys- We note that 
the CI velocities in gaps between filaments are generally 
redshifted compared to the velocities measured in the bright 
filaments themselves. We do not find evidence for the east- 
west galactic rotation seen in Hi emission. A careful com- 
parison of the particular regions covered by our SparsePak 



pointings to the H I velocity map of Miihle et al. ( 2005 their 
fig. 4), shows that in fact most of the equivalent Hi gas is 
either consistent with being at Vsya or very disturbed with 
no systematic pattern (see also Stil Sz Israel|2002 Heckman 
et al.|1995 L It is likely therefore that any signature of rota- 
tion in the ionized halo is either masked by the disturbed, 
strongly disrupted velocity field, or dominated by the effects 
of large-scale bubble expansion. 

The C2 radial velocity map (Fig. [t] right panel) shows 
the velocity variation of the second Gaussian component. As 
described in Section [2.2.11 this is by definition the fainter of 
the two brightest components fitted, and in all cases either 
represents a broad underlying component or the other half 
of a split-line profile, depending on what was required to 
fit the integrated line shape. Within the central ~500 pc, 
inside the broad component region, the velocities measured 
are very close to those found for CI. However at larger radii, 
where C2 is the fainter component of a split-line, we see ev- 
idence for ordered galactic-scale outflows. In the north, the 
blueshifted component is stronger (therefore in our conven- 
tion assigned to CI), and has velocities ranging between —20 
and kms~^ (relative to iisys), whereas the redshifted com- 
ponent (C2) has relative velocities of up to -1-150 kms~^. 
The maximum expansion velocity (~85 kms~^; calculated 
as half the difference between the radial velocities of the 
two individual components) is found to the north-west of 
the polar axis corresponding to the centre of shell B (see 
Fig. [2| ), and agrees well with the measurements of [Martin] 
(11998 1. In the north-east, a small number of split-line de- 
tections are coincident with the position of shell E (Fig. [2|, 
implying expansion velocities of up to 50 kms^^. In the 
south, the redshifted component is stronger (hence assigned 
to CI) and exhibits velocities between and -1-20 kms~^. 
The blueshifted component (C2) has corresponding radial 
velocities of up to —160 kms~^ (relative to «sys), giving a 



The galactic outflow in NGC 1569 9 





50 100 150 



I^^^^^^^^JBI^^B^^^^^^^^^^ 4.6 I I^^^^^^^B^^^^^^^BMI^^ 4.6 

Figure 5. H« flux maps for the four SparsePak fields for CI and C2 (left and right panels respectively). A scale bar is given for each 
plot (log units, arbitrary scale). Hatched spaxels represent non-detections, and the origin is the field centre of position 1(E). North is up 
and east is left. A reproduction of the WIYN Ha image is shown underneath the CI map for means of comparison. 
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Figure 6. FWHM maps of the Ha line for the four SparsePak fields, 
of kms~^, corrected for instrumental broadening. 



Left: CI and right: C2. A scale bar is given for each plot in units 



peak expansion velocity of 90 kms ^ to the south-west of 



SSC A, coincident with the centre of shell A (Fig. 



2 also 



agreeing well with the velocities quoted in Martin |1998 1 
Expansion velocities of >40 kms~^ are observed through- 
out the whole length of the cavity encompassed by shell A 
and the western arm. Blueshifted C2 velocities of —50 to 
—75 kms~^ are also seen in the northern half of position 



3(F), implying expansion velocities of up to 100 kms"'^ for 
shell F. We summarise the expansion velocity measurements 
for the four supershells in Table [2j together with their di- 
ameters and dynamical ages derived in Section [5. 2[ 



In three spaxels in position 2(B) (numbers 64, 73 and 
76), coinciding with the centre and western edge of shell 
B, we detect three distinct Ha line components similar to 
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Figure 7. Radial velocity maps of the Ha line for the four SparsePak fields. Left: CI and right: C2. A scale bar is given in units of 
kms~^, relative to fsys, where the vertical lines denote a velocity of 0. 
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Figure 8. Spatial distribution of electron density (in units of 1 



for the four SparsePak fields for left: CI and right: C2, calculated 



from the flux ratio of the [S ll] doublet. The colour map used is shown beneath each plot and ranges upwards from the low density limit 



of this indicator (~100 cm~^; |Osterbrock|l989t . 



the example shown in Fig. [4j;. The 'extra' component, with 
respect to the neighbouring double-component (narrow + 
broad) lines, lies at a velocity close to Ways, so could represent 
quiescent halo gas or be associated with part of an as yet 
undetected shell or bubble. 

The consistent reversal of strengths of the two major 
components of the split line from the north to south agrees 



well with previous studies ( [Heckman et al.]|1995[ |Martin| 
,1998i ), and demonstrates there is a definite preferred out- 
flow direction in both the northern and the southern hemi- 
spheres. This strengthens the hypothesis that the actual axis 
of the NGC 1569 galactic outflow may be inclined to the 



observer ( Stil & Israel 2002 I , in the sense that the north- 



ern outflow is tipped away from our sightline, rather than 
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simply that the bubbles are expanding out of the disc along 
random paths of least resistance. 

3.1 Electron Density 

Fig. |8] shows a map of the electron density, de- 
rived from the flux ratio of the [S ii] AA6717,6731 doublet, 
assuming ll] emission is detected in all 

spaxels that also have bright Ha emission, but often the 
densities implied wer e below the low d ensity limit for this 
indicator(~100 cm"''; Osterbrock| 1989 1 . Hence, we have set 



the minimum value of the colour range to 100 cm~ , mean- 
ing that spaxels where [S ll] emission was detected but im- 
ply densities below the low density limit are shown as black. 
Spajcels where [S ll] is not detected or where an unphysical 
flux ratio is measured (resulting easily from the large mea- 
surement errors) are represented by hatched lines. 

The majority of the ionized gas in the outer-wind re- 
gion is at or below the low density limit, as expected for a 
rarefied expanding outflow. Some regions in shell B [in po- 
sition 2(B)] and shell E (in position 1(E); see Figs [5] and |3| 
exhibit measurable densities in CI (left-hand panel) of or- 
der 500 cm~^. In parts of position 2(B), and a few spaxels 
in position 4(A), we are able to detect two line components 
and measure the electron densities of both to be of the or- 
der ~500 ± 300 cm~^. Unfortunately, the multi-component 
nature of the lines and the variable S/N have meant that 
the associated uncertainties on the measurements are high. 



3.2 Nebular diagnostics 

The flux ratios of [Sll]/Ha and [Nil] /Ha can be used as 
an indicators of the excitation level of the gas, and can be 
used to ascertain whether the gas is shock-excited, since the 
forbidden-line emission (particularly [Sll]) is enhanced in 



shocks (Dopita & Sutherland 1995 Calzetti et al. 20041 



Unfortunately, the strength of [Nil] in NGC 1569 is very 



weak compared to Ha (of order 2-5 per cent; see also Heck- 
man et aL||1995[ l due to the low metallicity of the gas, but 
nevertheless is detected in the brighter regions. 

Fig. [9] shows spatial maps of the 
[Sii](A6717-hA6731)/Ha and [Nii]A6583/Ha ratios. The 
limiting factor in both these maps is the S/N of the fainter 
forbidden-line emission; only where we have been able to 
detect emission have we been able to measure the line 
ratios. Having said that, the S/N of the forbidden lines is 
high enough in some places for multiple components to be 
detected, so to avoid confusion we use only the ratio of the 
brightest component (CI) identified by our Gaussian fitting 
routine for this ratio. 

The majority of spaxels where [S 11] emission is detected 
have a log([S ll]/Ha) « —0.5 to 0, and the highest ratios are 
generally associated with the faintest Ha emission whereas 
the lowest ratios are found where the Ha flux is highest. 
Spaxels with detectable [Nil] emission have a log([N ll]/Ha) 
~— 0.2 to— 1.3, and the lowest ratios are associated with the 
brightest Ha-emitting gas. As with [Sll]/Ha, this trend is 
correlated well with galactocentric distance, where the two 
ratios both increase from the region around SSC A towards 
the outer halo. 

These trends agree well with the measurements of |Heck- l 
man et al. (19951, who mapped the [Sll]/Ha and [Nil] /Ha 



ratios out to ~±50" from SSC A using two perpendicular 
slits, except that they found values of log([S ll]/Ha) « —1.2 
at the location of SSC A. However, [Heckman et al^ did not 
mention whether or not they fltted the [S II] lines with multi- 
ple components, so we do not know whether they measured 
the flux ratio of the integrated line-flux, or that of just the 
brightest component. If the former was true, then a signif- 
icant part of the discrepancy could arise from the varying 
contribution from the different line components. Further- 
more, most of the central region covered by their slits is not 
well sampled by our SparsePak flbres, and we may easily 
have missed the very low ratios they find in the inner 20" 
or so. 

A number of spaxels exhibit line ratios that are worth 
drawing attention to. Fig. [lO] shows two such spectra 
from the outer-wind regions with simultaneously very high 
[Sll]/Ha and very low [Nil] /Ha ratios. The upper spectrum 
has a [Sll]/Ha ratio of ~1.25 (log([S ll]/Ha) ~ 0.10), rep- 
resentative of the whole south-east of position 1(E) where 
most of the highest [Sll]/Ha ratios are found. Since there 
is no detectable [N 11] emission (see Fig. [9j right panel) , we 
can place an upper limit on the [Nll]/Ha ratio of <0.1. The 
lower spectrum shows a similar situation from the centre of 
position 3(F) (arc 1) where [Sll]/Ha is >0.5 and [Nll]/Ha 
is <0.1. High [Sll]/Ha ratios are also found clustered in the 
centre of position 2(B), position 3(F) (arc 1; lower spectrum 
in Fig. 10 1, and position 4(A). 



A ratio of log([S ll]/Ha) > lies in the shocked region of 
a [Olll]/H/3 vs. [Sll]/Ha diagram, for any sensible values of 
[Oiii]/H/3 dPopita fc Sutherland|1995l [Calzetti et al.|2004| ), 
thus indicating that at least shell E may be the result from 
an outward starburst-driven shock. Shocks are certainly ex- 
pected in the outfiow, since the hot superbubble will be ex- 
panding into preexisting Hi gas and/or other halo material, 
which is effectively at rest with respect to the outflow. The 
existence of shocks is also supported by the detection of 
strong X-ray emission from these shells (particularly shell E 
where the corresponding Ha emission is weak; [Martin et al.| 
Section 5.2 1. The shock excited component may be- 



2002 



come more obvious with increasing radius due to the fact 
that further from the central star clusters, the excitation of 
the gas is less dominated by the effects of photoionization. 



4 COMPARISON BETWEEN SPARSEPAK 
AND GEMINI GMOS/IFU DATA 

Before we discuss the implications of the broad component 
and what it may represent physically, we should compare its 
characteristics to the corresponding broad component found 
using our Gemini GMOS/IFU observations (Paper I Paper 

0- 

One way to compare the datasets quantitatively is by 
examining the flux ratio of C2/C1 - an important indica- 
tor of the signiflcance of the two component^ Although 
the quality of the data presented here is not high enough to 
identify any trends in this ratio across the observed flelds, 



^ We must be careful in diseussing this component since here 'C2' 
refers to both a broad underlying component and the other half 
of a split line profile, depending on context, whereas in Paper l| 
and |Paper II| 'C2' refers only to the broad underlying component. 
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Figure 9. Le/i; the spatial distribution of the [S ll](A6717+A6731)/Ha ratio and right: [N ll]A6583/Ha ratio for the four SparsePak fields 
(CI only). 
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Figure 10. Two examples of spectra that show high [Sll]/Ha 
ratios. The lower spectrum is from the centre of position 3(F) 
(fibre 25; log([S Il]/Ha) —0.25), and the upper spectrum (offset 
by 80 flux units) is from the far south-east of position 1(E) (fibre 
12; log([S Il]/H«) 0.20). Both are smoothed by 5 A, and shown 
on an arbitrary but relative flux scale. 



what does deserve attention is why we find such a high av- 
erage value of C2/C1 (~0.7) compared to the value found 
from the CMOS data (<0.3; [Pe^ern] ). 

It is simple to understand a high C2/C1 ratio for cases 
where C2 is one half of a split-line profile; in our assignment 
convention CI is by definition the brighter of the two compo- 
nents (C2/C1 < 1), and since in general the intensity of the 
two split-line components is comparable (to within a factor 
of a few), high C2/C1 ratios are expected. However, we also 



find C2/C1 ratios of ~0.2-0.7 within the broad component 
region, where from our CMOS data we find C2/C1 < 0.3. 

To investigate this discrepancy, we re-examined a num- 
ber of Ha profiles in the SparsePak spaxels closest to SSC 
A and the CMOS IFU fields (see inset to Fig.jsf, and found 
that in some of these cases a subtle extra component had 
been overlooked by the automatic fitting procedures. Fig. |ll| 
shows the Ha profile from fibre 50 of SparsePak position 
2(B) (the fibre nearest to CMOS IFU position 2 and cov- 
ering part of the bright Hll region No. 2; Waller 19911, 
and is the most obvious example of how, after close inspec- 
tion, a three-component fit is found to be more appropriate 
(shown in panel b) than the original, automatically deter- 
mined double-Gaussian fit (panel a). The residuals shown 
below each plot clearly demonstrate the improvement in 
fit quality, even though the formal improvement is not 
statistically significant. In this particular case, the effect of 
changing C2 of the double-Gaussian fit to C2 of the triple- 
Gaussian fit (now correctly assigned to a broad underlying 
component rather than a second narrow component) is to 
increase its FWHM by 50 kms~^ to ~150 kms^^, and to 
decrease the C2/C1 flux ratio from 0.69 to 0.17, now in 
much better agreement with that found from the GMOS 
data. This mis-fltting also offers an explanation for why the 
widths of C2 near SSC A measured from our SparsePak 
data are much less than the average valu es measu red from 
the GMOS observations (200-300 kms" 



Paper II I. 



Unfortunately, this kind of detailed fitting can only be 
done for the highest S/N spectra found near the centre of 
NGC 1569, so cannot be used reliably to track the C2/C1 
ratio over the full extent of the broad component region. 
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Figure 11. Ho profile from fibre 50 of SparsePak position 2(B) (nearest to GMOS IFU position 2 of |PaperT| see Fig.jsjl, illustrating how 
a more careful modelling of the line profile reveals that a three-component fit (panel b) characterises the integrated shape more accurately 
than the automatically determined double-Gaussian fit (panel a). The addition of a third component and the correct assignment of the 
profiles results in the FWHM of C2 increasing from 100 to 153 kms"'^ and the C2/C1 ratio changing from 0.69 to 0.17, now in much 
better agreement with the GMOS data for the same region ( [Paper Observed data is shown by a solid black line, individual Gaussian 
fits by dashed red lines (including the straight line continuum level fit), and the summed model profile in solid red. Flux units are 
arbitrary but on the same scale. The fit residuals are shown below the spectrum in both cases. 



5 DISCUSSION 

The central part of the ISM of NGC 1569 has been severely 
disturbed by the effects of an intense starburst event. 
Through deep Ha imaging and wide-field spatially resolved 
spectroscopy, we have investigated the state of the ionized 
gas in the halo of this galaxy. A detailed decomposition 
of the Ha line profiles shows that within a distinct region 
~700 X 500 pc in size and roughly centred on the location 
of SSC A, they are composed of a bright, narrow com- 
ponent with an underlying, broad component (FWHM < 
150 kms~^). We refer to this region as the broad component 
region. At larger radii, the broad component disappears and 
the profile shape becomes double-peaked, exhibiting charac- 
teristics of shell expansion. 



5.1 The inner galaxy and the broad component 



In IPaperT) we discussed in detail the possible mechanisms 
that could give rise to the underlying broad component 
found in the central regions of the galaxy. We concluded 
that the most likely cause is the impact of strong winds 
from the massive star clusters setting up turbulent mixing 
layers ( [Begelman &: Fabian] [T990l ) on the surfaces of ISM 
cloud clumps. Models predict that the ionized gas emission 
from these layers would exhibit line widths characteristic 
of the turbulent velocities within the layers ( [Slavin et al.| 
19931. The subsequent ablation and/or evaporation of ma- 
terial from these surfaces results in a pervasive, highly tur- 
bulent velocity field, explaining why the broad underlying 
emission is present throughout the inner regions of the bur- 
geoning outflow. 

Through the arguments presented in Section|4j we have 
established that the properties of the broad line component 
seen with SparsePak in the regions near to SSC A are equiv- 
alent to the broad line found in our Gemini GMOS/IFU ob- 



servations ( Paper I Paper II I . We can therefore confidently 



associate the two broad line components with each other and 
thus the physical mechanisms that produce them. However, 
the low spatial-resolution of the SparsePak array introduces 
a greater uncertainty to the degree unresolved kinematical 
components contribute to the line width. For example, one 
SparsePak fibre covers a similar area to the whole field of 
view (FoV) of the GMOS IFU (5 x 3.5 arcsecs; see Fig.|3|, 
and in |Paper II| we found C2 to have a radial velocity spread 
of >90 kms~ over just one IFU field. We must therefore 
accept that all our FWHM measurements will be artificially 
elevated by unresolved components at the level of a few tens 
of kms~^. 

An interesting and pertinent question arises from our 
study: why might the extent of the broad underlying emis- 
sion be limited to a specific area (the broad component re- 
gion), and what significance does this have? Our data show 
that the radius at which the broad line is no longer iden- 
tifiable also roughly corresponds to the point that the Ha 
profile starts showing a secondary narrow kinematic compo- 
nent. We speculate therefore that this radius marks a sig- 
nificant transition point in the development of the galactic 
outflow, where ordered expansion begins to dominate over 
turbulent motion, as is expected as a superbubble emerges 
from a disc. This possibly significant result warrants detailed 
further investigation at higher spatial-resolution and S /N to 
better quantify this transition region. 



5.2 The ionized halo 

The ionized gas morphology suggests a situation whereby 
individual bubbles, blown by a star cluster or SNe, expand 
into the halo through channels or 'chimneys' in the ISM 
distribution that were carved by previous generations of 
bubbles. Once a channel is cleared, swept up material can 
fiow through the same channel with much less energy losses, 
i.e. less mechanical energy is converted into thermal energy. 
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Figure 12. (a) Chandra medium-band X-ray contours ( [Martin et al.|[2002| l and (b) Hi column density contours | |Miilile et al.|[2"005[ l 
overlaid on the WIYN Ha image of Fig. [l] The shaded area on the bottom panel indicates the Hi 'hot-spot' (region of velocity crowding); 
a number of other features discussed in the text are also labelled. 
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In this interpretation, the outflow could be likened to a se- 
ries of cannon bursts, where each explosive event (possibly 
the formation of a massive star cluster or a collection of star 
clusters) has produced one of the bubbles we observe today. 

In Fig. |12[ we overplot contours representing medium- 
band X-ray emission (panel a) and H I column density (panel 
b) on the Ha image of Fig.[l] It is striking how well the X-ray 
emission morphology follows the Ha distribution (as noted 
by previous authors; [Heckman et al.1|1995| [Martin et al.| 
|2002[ ). However correspondences can also be found in the 
Hi map. 



Firstly, the Ha counterpart to the Hi chimney (Israel 
|fc van Driel||1990| |Muhle et aL||2005[ ) can be clearly seen 



extending to the south-east of SSC A (labelled on Fig. 12 a) 



This chimney appears to be the primary channel through 
which a significant amount of the disc ISM is escaping and 
forming superbubble F. Less prominent 'holes' in the Hi 
distribution also correspond well to the centres of shells A 
and E. Secondly, [Miihle et "aL| ( |2005[ ) identify a 'hot spot', 
a region of velocity crowding, on the western edge of the 
disc that they interpret as resulting from the impact of in- 
falling Hi gas from an interacting cloud. This area, shaded 
on Fig. |12^ , is just exterior to the peak of the Ha emission 



and the location of the large active Hll complex 2 (Waller 
|1991[ ). Multiwavelength observations of this part of the disc 



support the existence of a large molecular gas reservoir ( Tay 



lor et al. 1999 ) and the existence of continuing star- and 
cluster- formation ( [Tokura et al.|p006l [PapeFT] ), suggesting 
a possible reason for the irregular, off-centre location of the 
star burst and the two SSCs A and B. 

Perhaps more significantly, however, the hot-spot is lo- 
cated just to the west of filament 6 (the bright western arm, 
shell A), where a strong X-ray enhancement is also detected 
(Fig. [l2]a,; [Martin et"ar]|2002[ ). Why the western arm is so 
bright and part of such a large, coherent structure has been 
a standing problem since it was first detected ( Hodge|1974 |. 
We suggest that shell A was initially formed in a similar 
fashion to the other superbubbles, but as it expanded out 
into the halo it encountered a barrier formed by the infalling 
cold gas. The thermalisation of the outflowing material by 
the interaction resulted in increased radiative losses (show- 
ing up in both Ha and X-rays), hence the slowing of the 
bubble's expansion on the west side. The bubble then be- 
came distorted as the outflow was diverted from its original 
path, producing the structure that we see today. 

This leads us to consider the fate of the bubbles and the 



outflow itself. Martin ( 1998 1 interpreted the fact that the ve- 



locity separation of the two line components in shell A were 
still increasing when the fainter side faded below her detec- 
tion limit as evidence that this bubble may have ruptured. 
However, our deeper SparsePak observations show that this 
is not the case. We measure decreasing velocity separations 
in both shells A and B at the very edge of our combined 
FoV. In shell A, the expansion velocity peaks ~330 pc south 
of SSC A (at 90 kms"'^), then begins to decrease until we 
can no longer detect Ha emission. In shell B, the expan- 
sion velocity peaks ~375 pc north of SSC A (at 85 kms"'^), 
we again see a decrease in «cxp up to the edge of our FoV 
at distances of >1 kpc (although the S/N of our spectra 
at these radii is quite low). In support of this, signatures of 
shell blow-out are not seen in the Ha, H I ( [Miihle et al. 12005] ) 
or X-ray ( Martin et al.|200"2 | morphologies (Figs [Ij and 12 1. 



Firstly, there are no large-scale holes or under-densities in 
the H I column density of the halo that could be indicative 
of outflow channels. Secondly, the soft- and medium-band 
X-ray emission morphologies are very well matched to the 
Ha distribution, indicating that hot gas is not venting out 
of the bubble interiors (Fig. \V2^ \ [Martin et al.|2002 I. 

These conclusions do not, however, exclude the possi- 
bility of future blow-outs. The best estimates for th e escape 
speed in NGC 1569 are of the order of ~100 kms"^ (iMartin 



Fig.[2j 

ction p] 



19981, which when compared to our Ha expansion veloc- 
ity measurements, or to the kinetic temperature of the X- 
ray emitting gas, imply that the material in the shells may 
have the ability to escape the galaxy's gravitational poten- 
tial eventually. 

Measuring sizes of the four superbubbles from 
and using the expansion velocities calculated in Sectii 
we can estimate approximate dynamical ages for the shells. 
Shells A and B have radial extents of ~1.1 kpc and current 
expansion velocities on the order of 80-90 kms~^, giving 
dynamical ages of < 10-15 Myr. Estimating the true extents 
of superbubbles E and F is difficult considering their low 
surface brightnesses and more chaotic morphologies. Using 
the dashed lines plotted on Fig. [2] as a guide, we measure 
diameters for these bubbles of ~1.3 kpc, and assuming con- 
stant expansion velocities of ~50 kms~^, this gives ages of 
~25 Myr (or less if the shells were expanding faster in the 
past). These measurements are summarised in Table [5] 

We can compare these ages to the star and cluster for- 
mation histories derived in previous work. We know NGC 
1569 to have undergone a burst of ^constant rate star for- 
mation that stopped ~5-10 Myr ago from resolved stellar 
population modelling (Greggio et al. 19981. More specifi- 



cally however, we know from evolutionary synthesis mod- 
elling that SSCs A and B were formed ~ 10-30 Myr ago, 
and cluster 30 between ~30-100 Myr ago (Anders et al. 



|2004[ ). [Anders et al.| also concluded that the last major burst 
of cluster formation must have occurred some 25 Myr ago. 
These results fit in very well with our derived bubble ages, 
and lend support to the suggestion that each bubble forms 
after the creation of, for example, a massive star cluster. 



6 SUMMARY AND CONCLUSIONS 

We have presented WIYN MiniMo deep Ha imaging of NGC 
1569 covering a field-of-view of 9.5 x 10.7 arcmins. The depth 
and large dynamic range of the observations have enabled 
the identification of previously undetected faint ionized fil- 
aments in the halo and the study of the gas morphology 
right down into the bright, central regions of the disc. We 
have also presented WIYN SparsePak "formatted field unit" 
observations covering the outer galactic wind flow of NGC 
1569 in four pointings with integration times of 4-4.5 hours 
per field. The large diameter of the SparsePak fibres makes 
this instrument ideal for probing the faint ionized gas found 
in the halos of galaxies. This light-collecting power allowed 
us to choose a high-resolution spectrograph set-up, enabling 
us to characterise the line profile shapes of the important 
nebular diagnostic lines of Ha and [S ll] to an accuracy lim- 
ited only by the S/N achieved. We now summarise our main 
findings. 

• We find Ha emission out to radii of ~1.5 kpc from the 
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disc, and detect emission in almost every SparsePak fibre 
over the combined FoV. The presence of such an exten- 
sive system of ionized filaments results from the ongoing 
starburst that is supplying both mechanical energy to eject 
material and the Lyman continuum luminosity to keep it 
photoionized. 

• Through detailed Gaussian line fitting, we find that 
within a distinct region ~700 x 500 pc in size, roughly cen- 
tred on the location of SSC A, the nebular emission line pro- 
files are composed of a bright, narrow (FWHM < 70 kms~^) 
component with an underlying, broad component (FWHM 
~ 150 kms"'^). At larger radii, we find two narrow compo- 
nents to the Hq line, each representing one half of a split-line 
profile. 

• By comparing our results to observations of the central 
regions directly surrounding SSC A (Paper I Paper III, we 



conclude that the physical mechanisms that give rise to the 
underlying broad emission seen within this zone must be 
the same as within the regions directly surrounding SSC 
A sampled by Gemini GMOS/IFU observations. The broad 
emission is most likely to result from turbulent mixing layers 
on the surface of cool gas clumps set up by the impact of the 
hot, fast-flowing cluster winds, and from evaporation and/or 
ablation of material from the clumps. 

• The extent of this broad component region, coincident 
with the point at which we start observing signatures of 
large-scale bubble expansion, may indicate a transition point 
where ordered expansion begins to dominate over turbulent 
motion. Further observations are needed to investigate this 
in more detail. 

• By combining our deep Ha imaging and spectroscopy, 
we redefine the spatial extents of superbubbles A and B, and 



confirm their published ( Martin 1998 1 expansion velocities 
(90 and 85 kms~^, respectively). We estimate the dynamical 
ages of these bubbles to be ~10-15 Myr. Contrary to what 
has been previously suggested (Martin 1998), we find no 
kinematic or morphological evidence to suggest that either 
of these two superbubbles have ruptured and are venting 
their interiors into the galactic halo. 

• Our data indicate that the halo of NGC 1569 contains 
only 4 superbubbles. Following the terminology introduced 



by Martin ( 1998 1, we propose that her superbubble F should 
should encompass the whole south-eastern bubble complex, 



where the velocity ellipses identified by Martin ( 1998 1 and 
used to define shells G and F, are simply one level of a 'hi- 
erarchy of structure'. Furthermore, we propose that super- 
bubble E should encompass what were previously referred 
to as shells D and E and the large north-eastern X-ray spur 
( Martin et al.|2002| . 

• We derive new measurements of the expansion velocity, 
foxp (calculated from the difference in the radial velocities 
between the two split-line components), for the superbubble 
complexes E and F of iicxp ~ 50 kms~^ and <100 kms~^, 
respectively. Assuming a diameter of ~1.3 kpc for these two 
structures implies dynamical ages of <25 Myr. 

• The derived ages of the supershells are consistent with 



the recent cluster formation history of NGC 1569 (Anders 



et al. 20041, implying that each shell is associated with a 
specific star-forming event, such as a young massive star 
cluster that can provide a large mechanical luminosity from 
its many type H supernovae. 

• The consistent reversal of strengths between the blue 



and red components in the northern and southern outfiows 
provides evidence of preferred outflow directions approxi- 
mately perpendicular to the inclined and flattened Hi disc 
( StU fc Israel 2002) . 

• In addition to characterising the Ha line proflle, we 
have also measured [Sll] derived electron densities and 
[S ii]A6717-l-A6731/Ha line ratios. We flnd that much of the 
ionized gas in the broad component region and in the outer- 
wind regions is at or below the low density limit, as is ex- 
pected for a rarefied outflow. We find the highest [Sll]/Ha 
ratios are associated with the faintest Ha fluxes and the 
largest galactocentric distances. log([S ll]/Ha) ratios >0 are 
found in the superbubbles E, B and A, indicating that in 
these regions the gas emission may be signiflcantly shock- 
excited. 

In summary, the outflow in NGC 1569 appears to con- 
sist of several superbubbles in various phases of develop- 



ment, as noted by Martin (19981. This situation is further 



complicated by the disturbed state of the H I that includes 
features well out of the main plane of the galaxy (Fig. |12| 3). 
Our data confirm this model and thus indicate that the evo- 
lution of the outflow will be determined by the development 
of the superbubbles. In superbubbles A, B, and F ionized 
gas arcs seen in Ha and the X-ray morphology suggest that 
much of the hot gas still is confined, albeit moving at veloci- 
ties that are comparable to those needed to escape. Thus the 
NGC 1569 outfiow does not currently appear to be in the 
form of an approximately steady state galactic wind, even 
though it may eventually lead to mass loss from the system. 
It therefore differs from the well known M82 outflow, whose 
outer regions can be modelled by a supersonic galactic wind 



(e.g.|Suchkovetal.|1994| 


ShopbeU & Bland-Hawthorn|1998| 


Zirakashvili & V61k||2006 
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